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Introducción 


Una de las características más interesantes de la tele- 
detección es la posibilidad que proporciona el hecho 
de disponer de información periódica de un mismo 
lugar a lo largo del tiempo a un coste razonable. Esta 
propiedad facilita significativamente la realización de 
estudios multitemporales sobre las zonas de interés, 
con el fin de poner de relieve la evolución de ciertos 
rasgos en el territorio. De hecho, el análisis de cam- 
bios se ha convertido en una pieza fundamental den- 
tro de las aplicaciones de la teledetección. 


La identificación de los cambios entre dos imágenes 
de fechas diferentes ha sido abordada desde hace dé- 
cadas mediante el análisis visual. La fotointerpreta- 
ción tenía en ello uno de sus objetivos primordiales. 
El intérprete analizaba las formas, los colores, las tex- 
turas, los matices, la sombras y cuantos atributos vi- 
suales percibía y comparaba las dos imágenes para, 
finalmente, establecer los elementos que habían apa- 
recido, desaparecido o, simplemente, cambiado en la 
escena. 


El análisis digital de imagen ha venido a mejorar y 
profundizar en la identificación de cambios que in- 
cluso al operador podían escapársele. Sin embargo, 
aún hoy, carece del poder de síntesis que la mente 
humana desarrolla. Pequeñas diferencias entre las 
formas presentes en las imágenes o sutiles cambios 
de tonalidad pueden ser interpretados automática- 
mente como cambios, cuando un analista, en cambio, 
puede pasarlos por alto y considerar que no constitu- 
yen cambio en sí mismos. La comparación automática 
de imágenes no es inteligente hoy por hoy, por lo 
cual, las fuentes a comparar y necesitan cumplir unos 
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requisitos muy estrictos con carácter previo. 


Estos requisitos concretan en la invariabilidad de los 
elementos de imagen que realmente no hayan sufrido 
alteración con el paso del tiempo. Y dicha invariabili- 
dad debe ser satisfecha tanto desde el punto de vista 
gSeométrico como radiométrico. 


Como la comparación digital suele estar orientada bá- 
sicamente a la celda de imagen y, por lo general, el 
análisis automatizado no tiene en consideración fac- 
tores tales como el contexto, las formas, etc., en el 
análisis, pequeñas desviaciones de la posición de un 
mismo píxel serán interpretadas como un cambio. El 
ejemplo más paradigmático es el de una imagen sufi- 
cientemente abigarrada que no ha cambiado en abso- 
luto, pero que en una segunda fecha aparece despla- 
zada una sola celda en cualquier dirección con res- 
pecto a la original. El intérprete humano concluirá 
rápidamente que durante el plazo estudiado la es- 
cena no ha sufrido alteraciones, mientras que el aná- 
lisis digital advertirá de numerosos cambios. Por esta 
razón es de vital importancia la correcta y absoluta 
superponibilidad de las imágenes comparadas, lo 
cual requiere simplemente una eficaz corrección geo- 
métrica, incluida la ortorrectificación cuando sea ne- 
cesaria. 


De otra parte, se impone también un requisito adicio- 
nal, a saber, que la radiometría de las imágenes com- 
paradas sea equivalente, lo cual es necesario porque, 
de no serlo, cualquier desviación de los valores de 
luminancia, por ligera que fuera, podría ser interpre- 
tada como cambio. 
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Para identificar cambios temporales entre dos fechas 
con el objeto de actualizar una cartografía, por ejem- 
plo, o estudiar la evolución de un determinado fenó- 
meno, es necesario que los cambios reales en la ra- 
diancia medida de cada celda sean mayores que los 
debidos a factores externos como las diferencias at- 
mosféricas o el ángulo solar. El ejemplo extremo se- 
ría el de dos imágenes de un territorio que no ha 
cambiado en absoluto, pero que en una segunda fe- 
cha ha resultado afectado por una ligera calima ori- 
sinada por plumas de polvo en suspensión. Nueva- 
mente el sistema de tratamiento digital puede infor- 
mar de un cambio en toda la escena, aunque el intér- 
prete dirá que la segunda imagen es igual que la pri- 
mera, sólo que perturbada por el efecto atmosférico. 
En cualquier caso, para neutralizar las consecuencias 
de la variación de la irradiancia entre las fechas, lo 
que repercute en los valores de radiancia registrados, 
es preferible siempre realizar todos los análisis sobre 
reflectancia. 


Por todo lo expuesto, es obligada la corrección radio- 
métrica rigurosa basada en modelos como 
LOWTRAN, MODTRAN, FLAASH, etc. Para imáge- 


Métodos de normalización 


Como se ha dicho, es un requisito previo que las imá- 
genes a comparar sean equivalentes desde el punto 
de vista radiométrico, de tal manera que las variacio- 
nes en los ND de las celdas entre una y otra fecha, 
más allá de la variación inherente a la propia cu- 
bierta, sea debido a un cambio real en la ocupación 
de suelo. Aun habiendo realizado una normalización 
topográfica y la corrección atmosférica previa, suelen 
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nes aéreas, con menor incidencia del efecto atmosfé- 
rico, en muchos casos es suficiente realizar una nor- 
malización radiométrica lineal, alterando la reflectan- 
cia de una de las dos imágenes reescalándola y su- 
mando o restando alguna constante, de modo que la 
imagen de referencia y la normalizada sean equiva- 
lente desde el punto de vista radiométrico: 


P =44Pi +D, (1) 


donde p, es la reflectancia en la banda y p, es el 


valor normalizado de reflectancia normalizada que se 
pretende conseguir en esa banda. De este modo se 
altera la imagen objetivo con la ayuda de los paráme- 
tros ax y bx, que se extraen de la imagen de referen- 
cia. 


Dado que la invariabilidad geométrica se confía a la 
corrección geométrica, no se insistirá más en ello en 
este capítulo. Quedan por establecer los procedimien- 
tos de homologación de las imágenes desde el punto 
de vista radiométrico, a lo cual se dedica el siguiente 
apartado, para pasar finalmente a exponer algunas de 
las técnicas más usuales en la detección de cambios. 


radiométrica 


persistir diferencias ligeras en los ND de celdas su- 
puestamente invariantes. Es por esto por lo que se 
impone la necesidad de normalizar una de las imáge- 
nes a los niveles radiométricos de la otra, la cual será 
considerada como referencia, siendo irrelevante cuál 
de las dos sea tomada como tal. En la exposición que 
sigue, se toma la imagen de la primera fecha como 
referencia, pero bien podría ser al contrario. 





Esquema del proceso de normalización 
radiométrica. 


El procedimiento general a seguir consistirá en una 
transformación lineal, calculada banda a banda. Para 
una determinada banda k: 


y=ay+b (2) 


en cuya expresión se considera la siguiente termino- 
logía, que será la que se mantendrá en lo sucesivo: 


— ay b son dos parámetros deducidos de ciertos va- 
lores de la variable x, que representa el ND, o la 
reflectancia en su caso, en la fecha 1, que se toma 
como imagen de referencia. 

— y esel ND, o la reflectancia en su caso, en la fecha 
2. Corresponde a la imagen objetivo. 

— yes el ND normalizado, o la reflectancia normali- 


zada en su caso, en la fecha 2, con los valores de 
referencia de la fecha 1. 


Los distintos procedimientos de normalización sólo 
difieren en la forma de calcular a y b. a partir de los 
valores de la imagen de referencia x, que suelen de- 
ducirse de los estadísticos de esta imagen. 


Normalización por regresión en zonas 
invariantes 


Este método, conocido como No-Change Regression 
Normalization, se basa en la búsqueda de regiones 
de no-cambio de tierra firme y agua en el diagrama 
de dispersión multitemporal del IRp. 


Los parámetros a y b se calculan del siguiente modo: 





(NC 
AE 
(NC 
O yy 
¡q ad == a 
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de tal modo que la transformación es: 





ID Ar (3) 


Nótese que esta expresión tiene la forma de (1). En 
todas las ecuaciones anteriores: 


ez Lo 3 (x o (y yo ] es la covarianza 


entre las dos fechas para las zonas invariantes. 


l _ (NO Y? 
E S a yy ne] =0 “es la varianza de 
NC NC 


las zonas invariantes. 


Normalización por regresión simple 


Conocido como Simple Regression Normalization, 
éste es un método similar al anterior, pero en el que 
se extraen los estadísticos de toda la imagen, en lugar 
de hacerlo sobre una zona concreta de no-cambio. 





Ahora: 
o] 

a== 
O, 

b=x-ay 

y 

A 6 == 

Y=— (yy) +x (4) 
O y 

donde 
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Xy 


O. = Pl) (xy) es la covarianza [fecha 1- 
N 


fecha 2] de las dos imágenes. 


l 
O, = A — y) es la varianza de la imagen de la 
N 


fecha objetivo (fecha 2). 


Normalización por conjunto de control 
radiométrico 


La Radiometric Control Set Normalization utiliza la 
media aritmética de un subconjunto de celdas oscu- 
ras y brillantes extraído de las imágenes referencia y 
objetivo mediante la transformación tasseled cap. Es- 
tablecidos ambos conjuntos de datos, se calculan los 
parámetros con el siguiente procedimiento: 


qe 0 
(dl = 
751 
pe a ye e ye 
E 
y —y 


donde el superíndice (B indica las regiones de alto 
brillo (bright) y D indica las oscuras (dark). La nor- 
malización se lleva a cabo con: 


Py yO) (y+ y") 

LAS A.  —* (5) 
y Ty 

Normalización por corrección de la neblina 


En realidad, la Haze Correction Normalization, 
como se conoce en los textos ingleses, se refiere a 
cualquier tipo de efecto atmosférico dispersivo, bien 
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causado por aerosoles líquidos, como las nieblas, 
bien sólidos, como las calimas. Se basa en el siguiente 
postulado: sobre un área de reflectancia nula, la ima- 
sen debería aparecer oscura. De apreciarse alguna 
reflectancia aparente sobre este tipo de zonas, con de 
ND significativamente distintos de cero, debe ser la 
atmósfera la que está proporcionando ese nivel de 
brillo como consecuencia del efecto atmosférico, par- 
ticularmente debido a fenómenos de dispersión aero- 
sólica. 


Esta corrección se basa en hacer el mínimo ND de la 
imagen objetivo igual al de la de referencia. Los pa- 
rámetros son triviales: 


a=1l 


b sl A mín o má 
quedando: 
y Ñ y o Vi + A mín (6) 


Solamente desplaza el histograma de la imagen obje- 
tivo, sin escalarlo, para hacer coincidir los valores mí- 
nimos de los dos histogramas. 


Normalización por la media y desviación 
típica 

Conocido como Mean-Standard Deviation Normali- 
zacion, mediante este procedimiento se igualan las 
medias y las desviaciones estándar de la imagen ob- 
jetivo a las de la imagen de referencia. La ganancia y 
el sesgo de la transformación son: 





a= 
0 
b=x-ay 





Y==(y-y)+x (7) 


Se desplaza y reescala el histograma de la imagen 
objetivo para asemejarlo al de la de referencia. En 
éste, como en todos los procedimientos similares que 
pretenden igualar (matching) dos histogramas, se si- 
Sue el esquema siguiente: 


1. Se desplaza el histograma de la imagen obje- 
tivo hasta centrar su media en el origen 


(y-y). 

2. Se dividen los valores por la varianza o la des- 
viación típica del histograma objetivo y se 
multiplica por la del histograma de referencia 
que se desea obtener o, /0, . 


3. Se desplaza el histograma resultante hasta al- 
canzar su media el valor promedio de la ima- 
sen de referencia (+x ). 


Normalización por mínimo y máximo 


La Minimum-Maximum  Normalization ajusta, 
banda a banda, el contraste de la imagen objetivo 
para conseguir histogramas con los mismos valores 
máximos y mínimos de las bandas homólogas de la 
imagen de referencia. 


Métodos de normalización radiométrica 


A náx Xy ín 


más a mía 
b A mín O mín 


(al = 


La imagen normalizada se obtiene aplicando: 


A máx o A nín ( 


y pe Y mín ) a A mín (8) 
Vs o DÍ 


ve 


La transformación escala el histograma de la segunda 
fecha utilizando el criterio de valores mínimo y má- 
ximo, sin reparar en la varianza. 


Normalización por regiones 
pseudoinvariantes 


En esta transformación, conocida como Pseudo-Inva- 
riant Normalization, se analizan elementos de reflec- 
tancia invariante (Pseudo- Invariant Features, PIFSs). 
Su desarrollo es similar al visto en (7), pero en esta 
ocasión, los estadísticos se extraen solamente de un 
subconjunto de celdas de las que se tiene certeza de 
que su reflectancia no cambia con el tiempo, ni por 
la estación del año, ni por la hora, ni por el grado de 
humedad. Los parámetros de la transformación se ba- 
san en el cálculo de la media y varianza de esas zo- 
nas, igualando las de la imagen objetivo a las de la 
de referencia. 








O, 
a = e 
O, 
b=x" gy" 
A ¡0 dl —( pi —Í pi 
>= (y 3 )+z (9) 
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donde el superíndice (pi hace referencia a las regio- 
nes pseudoinvariantes. Para la aplicación del método 
pueden utilizarse superficies como el hormigón, las 


superficies de aglomerado asfáltico o ciertas cubier- 
tas de edificios, en general, elementos exentos de ve- 
setación y de materiales que cambien su contenido 
de humedad. 


Técnicas de detección de cambios 


En este apartado se recogen algunos de los procedi- 
mientos habituales para comparar dos imágenes ad- 
quiridas en fechas diferentes. El análisis multitempo- 
ral puede exigir un estudio más minucioso que reco- 
rra varias etapas entre el instante inicial y el final. De 
hecho, la evolución de una superficie puede requerir 
muchas observaciones a lo largo de un plazo amplio. 
Esto entraña cierta dificultad a la hora de representar 
la evolución, salvo que se emplee el tiempo en la pro- 
pia representación bajo alguna modalidad de anima- 
ción o se represente en forma de gráfica de dos ejes 
la evolución de alguna característica de una sola 
celda o de un conjunto de ellas, usando la variable 
tiempo como abscisa. Las explicaciones siguientes 
son válidas para la comparación de sólo dos momen- 
tos, que pueden ser el inicial y el final o bien cual. 
quier otro intervalo temporal intermedio. 


Composiciones en color multitemporal 


Cuando se desea comparar mediante un análisis vi- 
sual rápido la variación de una superficie entre dos 
fechas se puede recurrir a la composición en color 
RGB de una banda seleccionada de la imagen en cada 
una de las fechas. Como en este caso se dispone de 
tres canales, pero sólo de dos fuentes de información, 
una de las bandas ha de repetirse en la composición. 
Así pues, se selecciona la banda en la cual la reflec- 
tancia sea más variable entre las dos fechas. 
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Por ejemplo, si se trata de identificar la presencia de 
un cierto cultivo y se desea conocer cuándo ha estado 
presente sobre el terreno, se puede seleccionar la 
banda del infrarrojo próximo, en la cual el vigor ve- 
setal proporcionará crecientes niveles de reflectan- 
cia. Si se asigna la banda de la fecha inicial al color 
rojo (R) de la composición en falso color y la de la 
fecha final tanto al verde (G) como al azul (B), un 
incremento de la reflectancia entre ambas fechas 


Banda k — fecha 2 


AA A 


A A 





Composición RGB 





No cambio 


Disminución de la reflectancia 


Aumento de la reflectancia 


Composición en color multitemporal utilizando sólo dos fechas. 


Banda lRp — fecha 3 


Banda IRp — fecha 2 


Banda |Rp — fecha 1 





Composición RGB 





mostrará tonos en los que estará ausente la compo- 
nente R y predominarán las componentes G y B. El 
resultado, pues, serán tonos próximos al cian (C). 
Esto viene a indicar que, si el cultivo buscado se desa- 
rrolla entre las dos fechas y aún está vigoroso du- 
rante la segunda, sólo las celdas de color cian serán 
las que cumplan esa condición. 


Puede suceder que las parcelas con el cultivo bus- 
cado estuvieran en pleno desarrollo vegetativo en la 


primera fecha, pero en la segunda el cultivo o bien 
esté agostado o directamente haya sido cosechado y 


El fenómeno estudiado: 


No ha estado nunca 





Desaparece antes de la fecha 2 


Aparece tras la fecha 1 y 
desaparece antes de la fecha 3 


Aparece tras la fecha 2 


Está presente desde el inicio y 
desaparece antes de la fecha 3 


Aparece tras la fecha 1 y no 
desaparece posteriormente 


Está siempre 


A 


Composición en color multitemporal utilizando tres fechas. 
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labrado. En ese caso predominará el color rojo sobre 
la mezcla de verde y azul, aunque esta última no des- 
aparezca del todo. El color resultante será un rojizo 
más o menos puro. 


El tercer caso que puede darse es cuando el cultivo 
haya estado presente en las dos fechas consideradas, 
y con un vigor vegetal similar, o bien no haya estado 
nunca presente. En ambos casos los niveles asigna- 
dos a los canales R, G y B será similares y el resultado 
es que las parcelas que cumplan esa condición exhi- 
birán una tonalidad gris-parda. 


Sin embargo, dado que se dispone de tres canales, es 
posible introducir en el estudio multitemporal hasta 
tres fechas sucesivas y sintetizar el resultado en una 
sola imagen. El intérprete adiestrado identificará sin 
problemas la presencia, ausencia o evolución de un 
cierto cultivo, dentro de un conjunto reducido y co- 
nocido de ellos. En ese caso se asigna la banda del 
infrarrojo próximo de la fecha inicial, 1, al canal R, la 
de la fecha intermedia, 2, al G y la de la fecha final, 
3, al B. 


Recordando la teoría de la síntesis aditiva de forma- 
ción del color en el sistema RGB, la interpretación de 
los posibles casos viene dada así: 


— En las celdas con tonalidades gris o pardo oscuro 
el cultivo buscado no ha estado presente nunca. 
El suelo está desnudo. 

— En las celdas rojizas, el cultivo está presente en la 
fecha 1 pero desaparece antes de la fecha 2. 

— En las celdas verdosas, el cultivo aparece tras la 
fecha 1 y desaparece antes de la fecha 3. 

— El color azulado indica que el cultivo sólo aparece 
tras la fecha 2, y sigue presente en la fecha 3. 

— Si el color de la celda es amarillento, el cultivo 
está presente desde el inicio y durante la fecha 2, 
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pero desaparece antes de la fecha 3. 

= Si el color tiende al cian, el cultivo no está pre- 
sente en el inicio del período estudiado, sino que 
nace tras esa fecha, pero ya está desarrollado en 
la fecha 2 y se mantiene sobre el terreno en la 
última fecha. 

— Finalmente, si el tono es muy claro, gris o pardo 
brillante, el cultivo ha estado sobre el terreno du- 
rante todo el intervalo estudiado. 


Esta metodología se utiliza habitualmente en el se- 
guimiento de la ocupación de suelo por determinados 
cultivos beneficiarios de subvenciones por las admi- 
nistraciones públicas. En su origen se realizaba un 
análisis visual y hoy día se lleva a cabo mediante aná- 
lisis automatizado digital de imágenes. 


Diferencia entre imágenes 


Una de las opciones más inmediatas para detectar 
cualquier cambio entre dos fechas en una escena es 
restar las imágenes de reflectancia de ambas: 


DI=p, -Pp, +C 


Es obvio que, admitiendo un mínimo grado de varia- 
ción de la variable aun no habiendo cambiado la co- 
bertura del suelo, tras la diferencia de las imágenes, 
las celdas que no hayan mutado ostentarán unos va- 
lores próximos a cero. Por el contrario, las áreas de 
cambio mostrarán resultados sensiblemente positivos 
o negativos. 


Habitualmente, será necesario reescalar el contraste 
de la imagen diferencia, o trasladar su histograma, 
para hacerla apreciable en un sistema de visualiza- 
ción de 8 bits o superior. Ese histograma de los valo- 
res de la imagen diferencia tenderá a distribuirse 
sSaussianamente. Para definir las zonas de cambio en 
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- 255 255 


la imagen se deberá establecer un valor umbral para 
la frontera cambio-no cambio en el histograma aten- 
diendo, si es posible, a algún rasgo peculiar evidente, 
a izquierda y derecha del valor modal. La elección de 
este umbral es determinante para obtener resultados 
demasiado restrictivos si el valor está próximo a la 
moda de la distribución, esto es, perdiendo cambios 
y sólo identificando los muy ostensibles, con gran 





Histograma fecha 1 






XT 


Histograma diferencia 


Area de no cambio 


Histograma fecha 2 


Histograma reescalado 


Identificación de las áreas de no cambio mediante la diferencia de imágenes. 








255 


error de omisión, o bien obtener resultados muy la- 
xos, resultando cambios en celdas cuya reflectancia 
ha variado muy ligeramente, posiblemente como con- 
secuencia de la propia variación interna de la clase 
que representa y cometiendo así numerosos errores 
de comisión. 


Cocientes multitemporales 


Otra opción sencilla es dividir celda a celda los valo- 
res de reflectancia de una determinada banda de las 
imágenes de las dos fechas. 


cm=*2.0 
P; 


El valor de CM se mantendrá próximo al valor de C 
para las zonas invariantes. Este procedimiento tiene 
la ventaja de reflejar el valor relativo del cambio, 
frente al valor absoluto que proporciona el análisis 
por diferencias de imágenes. 


Para evitar valores fuera de rango y hacer compara- 
bles los cocientes, puede emplearse el índice del co- 
ciente multitemporal de diferencia normalizada, que 
aporta la cualidad de estar acotado entre 0 y 1. Para 
un sistema de visualización de 8 bits de profundidad: 


NDRI=| P2Pr 
Pp, FP, 
La interpretación del NDRI es la siguiente: 


— Si p crece mucho, el VDRI tiende a 1. 
— Si p, varía muy poco el NDRI está próximo a 0. 
—= Si p disminuye mucho, el NDRI tiende a —1. 
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También puede utilizarse una versión escalada del 
mismo índice: 


Ahora: 


— Si p crece mucho, el VDRT' tiende a 255. 
— Si p varía muy poco, el NDRT' oscila sobre 127. 
= Si p disminuye mucho, el NDRT' tiende a 0. 


Diferencia entre los índices de vegetación 


Al igual que se hace con las imágenes, se puede cal. 
cular el cambio que el índice de vegetación de dife- 
rencia normalizada ha sufrido entre las dos fechas 
sustrayendo uno de otro. 


DIV =NDVI, — NDVI, +C 


Puede usarse cualquier otro índice de vegetación, 
por ejemplo, el SAVI: 


DIV! =SAVI, —SAVI, +C 


De no añadir la constante C, ambos índices estarán 
acotados entre -2 y 2. Al igual que en la diferencia 
de imágenes, los valores de la diferencia de índices 
suelen presentar una distribución normal. Las zonas 
de no cambio tendrán valores próximos a 0 y las que 
hayan cambiado serán las de mayor valor absoluto, 
tanto positivo como negativo. Es conveniente buscar 
algún rasgo significativo en el histograma de la dis- 
tribución para delimitar el umbral de valores de la 
zona de cambio. 
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También se puede recurrir a la diferencia normali- 
zada de los índices de ambas fechas, obteniéndose el 
índice multitemporal de diferencia normalizada: 


NDVI, — NDVI, 


MNDVI = == ——_4 
NDVI, + NDVI, 


que proporciona valores normalizados entre —1 y 1. 


Otra forma de enfocar el análisis multitemporal con 
los índices de vegetación es aplicando el color tal 
como se ha explicado arriba sintetizando una sola 
imagen analítica con la información del índice en 
hasta tres fechas diferentes. Como antes, se puede 
utilizar el NDVI o bien el SAVI, si la escena muestra 
gran proporción de suelo desnudo. Para la composi- 
ción en color se asigna el cañón R al índice de la 
primera fecha, el G al de la segunda y el B al de la 
última fecha. La composición RGB resultante permite 
identificar las masas vegetales que aparecieron o des- 
aparecieron a lo largo plazo considerado, conside- 
rando la interpretación de los colores descrita ante- 
riormente. 


Regresión 


Se considera que en un contexto de no-cambio, la re- 
flectancia de la celda en el instante f2 es una función 
lineal de la correspondiente al instante tf para cada 
banda k: 


P, (1,)=a-p, (1,)+b 
siendo p, (t,) la reflectancia estimada en £2 a partir 
de la del momento fi y a y b, los parámetros del 


ajuste. El cambio real sufrido entre ambas fechas se 
calcula como la diferencia entre el valor realmente 
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observado y el estimado por regresión, que equival- 
dría al no-cambio: 


Análisis del vector de cambio 


Cada celda de una imagen multiespectral de k bandas 
está definida por su vector de características en un 
espacio vectorial k-dimensional. En el instante ft, el 
vector de características habrá variado con respecto 
al instante fi. El cambio se puede estimar midiendo 
la diferencia entre los vectores de características o el 
ángulo espectral entre ambos. 


Atendiendo al primer criterio, la magnitud del cambio 
viene determinada por la distancia euclídea espectral 
entre los extremos de los dos vectores de caracterís- 
ticas: 








midiendo el módulo del vector diferencia entre las 
dos posiciones espectrales de la celda en cuestión. 


Ángulo espectral de cambio 


Una forma alternativa de poner de manifiesto el cam- 
bio de la posición espectral de una celda es medir el 
ángulo formado por el vector de características de la 
celda entre los instantes inicial y final. 


La abundancia de bandas asegura no detectar falsos 
negativos, ya que aunque el ángulo espectral sea pró- 
ximo a O en algunas bandas, casi siempre se encon- 
trará alguna otra en la que esto no se cumpla y, en 


Momento 
La masa vegetal: 


Desaparece antes de t, 
Desaparece antes de t, 
Está siempre presente 
Aparece tras t, y desaparece antes de t, 
Aparece después de t, 


Aparece depués de t, 


— 
———— 
———_— 
— 
=> 
— 


|| E [4 5] mM mí Color predominante 


Munca ha estado presente 


Interpretación del color en la composición 
de índices multitemporales. 


ND 


Banda L / da , ¡e 
j h 





HD 
Banda K 


Diferencia entre los vectores 
de características en dos momentos. 


banda | 


Na 














Ángulo espectral de cambio. 


consecuencia, el ángulo espectral final será diferente 
de 0, indicando ello que habrá habido algún cambio. 
Solamente no se considerará cambio cuando el án- 
sSulo espectral continúe siendo O o muy próximo a 
este valor, tomado el conjunto completo de bandas, 
En este caso la celda habrá podido migrar de posición 
en el espacio espectral pero sólo en la dirección del 
origen de coordenadas espectrales. Lo cual indica 
una pérdida generalizada de la reflectancia, si es ésta 
la variable medida, que suele venir asociada a un au- 
mento de la humedad. Si lo que se registra es una 
pérdida de radiancia, puede estar relacionada tam- 
bién con una peor iluminación de la ladera. 


El ángulo espectral viene establecido: 


QA = arccos 





donde z, y z, son las componentes del vector de 
fl t) 


características de la celda en el momento 1 y 2, res- 
pectivamente. 


Transformación a componentes principales 


La información contenida en una imagen multiespec- 
tral está distribuida entre el conjunto de bandas, si 
bien no todas ellas proporcionan información inde- 
pendiente. Por el contrario, la información contenida 
en cualquier banda suele ser redundante en cierta 
medida con la contenida en otras bandas, sobre todo 
en las más cercanas en el espectro electromagnético. 
La redundancia en la información se manifiesta en 
una alta correlación entre los niveles digitales de cier- 
tos canales. 


Técnicas de detección de cambios MP 


La transformación o análisis de componentes princi- 
pales (ACP) proporciona unas nuevas bandas lineal. 
mente independientes mediante la combinación li- 
neal de las originales. 


La aplicación de esta metodología al análisis multi- 
temporal consiste en realizar la transformación con 
los dos juegos completos de bandas de las fechas a 
comparar. De esta manera, se extraen componentes 
principales, ordenados por varianza decreciente, del 
conjunto de bandas de las dos fechas. 


La mayor varianza corresponde a los primeros com- 
ponentes principales, que son los que contienen ma- 
yor cantidad de información. Han de contener, pues, 
aquellos rasgos presentes en ambas fechas, puesto 
que explican la mayor parte de información del do- 
blete. En ellos habrá que buscar los elementos inva- 
riantes. Los componentes principales de ordinal su- 
perior muestran las superficies que hayan sufrido 
cambios. Finalmente, los últimos componentes prin- 
cipales suelen retener solamente el ruido de la ima- 
gen y no representan necesariamente los cambios ne- 
tos reales. 


Una técnica alternativa que utiliza el ACP consiste en 
restar los componentes principales extraídos separa- 
damente de cada fecha. 


Correlación cruzada 


Se basa en la comparación de una imagen multies- 
pectral del momento f2 con una imagen clasificada en 
el momento f4. 


Se calcula el cuadrado de la diferencia entre el nivel 


digital de una celda en el instante £2 y el promedio de 
todas las celdas que, en ese mismo instante, están en- 
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globadas bajo el perímetro de clase a la cual pertene- 
cía en f1, según una clasificación realizada en esa pri- 
mera fecha. Se divide el valor anterior por la varianza 
de cada área y se suman los valores obtenidos para 
todas las bandas: 


Valores altos de Zc indican mayor probabilidad de 
cambio entre la clase inicial en t1 y la posterior en ft. 


Comparación posclasificación 


Este procedimiento es uno de los que goza de mayor 
difusión, en parte por la robustez de sus resultados y 
en parte por la facilidad en su interpretación. 


Como en todos los métodos, ha de comenzar con la 
seorreferenciación rigurosa de las imágenes de las 
dos fechas. A partir de entonces, se clasifican ambas 
con el mismo algoritmo. Es conveniente que la le- 
yenda sea idéntica, para lo cual, si es necesario, 
puede incorporarse la categoría no clasificado. 


Tras la clasificación, se comparan celda a celda las 
dos imágenes y se confecciona una tabla de contin- 
sencia con los resultados. Pueden ponerse en colum- 
nas las clases en el momento t1 y en filas las del mo- 
mento f2, de tal forma que cada celda de la matriz 
indica el número de celdas de la imagen que en la 
fecha inicial pertenecían a la categoría señalada en el 
encabezado de su columna y en la fecha final lo están 
en la categoría marcada por su fila. Como es obvio, 
no habrán cambiado de categoría las cantidades de 
celdas recogidas en la diagonal principal de la matriz, 
cada una en la clase que le corresponda por columna 
o línea. La suma de las celdas fuera de la diagonal 
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principal indica el número de celdas de imagen que 
han sufrido algún tipo de cambio entre las fechas exa- 
minadas. 


El estudio analítico de los cambios es sencillo. Sin 
embargo, es más compleja su presentación en forma 
de una imagen sinóptica que dé una idea fácilmente 
comprensible de las zonas de la escena que han per- 
manecido inalteradas y aquellas otras que han su- 
frido cambio, y de qué naturaleza ha sido éste. Hay 
que tener en cuenta que, si el número de clases con- 
siderado es n, si bien algunas mutaciones en la ma- 
yoría de los casos son inusuales o incluso inviables, 
el posible número de cambios podría llegar a ser n?. 
Una manera rápida de confeccionar una leyenda de 
n* categorías sin confusión entre todas las posibles 
mutaciones de clase y que mantiene la identidad de 
las celdas que no han cambiado es la siguiente: 


Si n2 es el ordinal de la clase en £ y ni el co- 
rrespondiente al momento fi, una asignación 
inequívoca podría consistir en etiquetar el 
cambio ni—>n2 calculando: 


(a, +1) +I=m | 


Obtenida la etiqueta del cambio, sólo resta 
asignarle un color de la paleta cualitativa pre- 
viamente confeccionada. 


No obstante, es discutible que la imagen sintética de 
cambios obtenida pueda ser fácilmente interpretable 
sin recurrir continuamente a consultar la leyenda. 
Por encima de 4 clases no es práctico representar los 
cambios en un solo mapa cualitativo, ya que la lec- 
tura de los colores o de las tramas de la paleta puede 
resultar difícil. Por ello, es más práctico, aunque 
requiera presentar un mayor número de imágenes, 
realizar n mapas, uno por clase, siendo n el número 





Imagen de 


1200x800 Clase1 Clase 2 Clase 3 Clase 4 Clase 5 


Clase 1 4 3 2 1 0 
Clase 2 9 8 7 6 5 
Clase 3 16 15 14 13 12 
Clase 4 25 24 23 22 21 
Clase 5 36 35 34 33 32 


Codificación de cambios para una imagen con 5 cla- 
ses iniciales y finales. Se confeccionan primera- 
mente los ordinales modificados de las filas (en 

blanco) a partir de los ordinarios de clase mediante 
(n2+1)? y se le sustrae el ordinal de la columna 

(también en blanco). Con ello se obtienen índices bi- 

unívocos de cambio de clase para las 25 posibilida- 

des, con los que confeccionar una leyenda. 


de clases, cada uno de los cuales contiene un máximo 
de n categorías, en los cuales se presenta la evolución 
de cada una de las clases de la leyenda. Puede 
realizarse desde la perspectiva de las clases 
originales, es decir, la superficie que en f¡ era clase c; 
ha pasado a ser c1, Ca, etc. (de ahí que se originen 
hasta n categorías), o también desde la perspectiva 
de la ocupación actual del suelo, o sea, lo que en t> 
es realmente c;, en ti fue c1, Ca, etc. Las aplicaciones 
informáticas de uso común en  teledetección 
proporcionan una imagen cualitativa multiespectral 
con n bandas, una por clase final, cada una de las 
cuales tiene a su vez n categorías, una por clase 
inicial, y son susceptibles de invertir de igual forma 
la dirección temporal del análsis. 


Clasificación multitemporal 


Este procedimiento clasifica como un solo conjunto 
los dos juegos de bandas correspondientes a ambas 
fechas. En el espacio de características se formarán 
dos tipos de nubes de puntos: aquellos que corres- 
ponden a celdas invariantes, que verán duplicados 
sus vectores de características, al haberse mantenido 
estáticas en el espacio espectral, y aquellas otras que, 
al haber sufrido cambios, presentarán un desdobla- 
miento de tales vectores. En el proceso de clasifica- 
ción no supervisada, será fácil identificar esos clus- 
ters mutantes. 


Técnicas de detección de cambios ML 
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